A1 = readmatrix('NiI2_freestanding_20nm_IV_30K.txt')
A2 = readmatrix('NiI2_freestanding_20nm_IV_40K.txt')
A3 = readmatrix('NiI2_freestanding_20nm_IV_50K.txt')
A4 = readmatrix('NiI2_freestanding_20nm_IV_60K.txt')

x1 = A1(:,1)
y1 = A1(:,2)

x2 = A2(:,1)
y2 = A2(:,2)

x3 = A3(:,1)
y3 = A3(:,2)

x4 = A4(:,1)
y4 = A4(:,2)

subplot(4,1,4,'Position',[0.15 0.1 0.8 0.2])
hold on
scatter(x1,y1,3,'filled','MarkerFaceColor',[0.2 0 0.8])
plot(x1,y1,'Linewidth',1,'Color',[0.2 0 0.8])
xlim([-13,13])
ylim([-3,3])
xticks(-12:4:12)
yticks(-2:1:2)
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
xlabel('Electric Field [MV/m]')
dummyh = line(nan, nan, 'Linestyle', 'none', 'Marker', 'none', 'Color', 'none')
%set(gca,'XMinorTick','on','YMinorTick','on')
legend(dummyh, '30 K')
legend boxoff
hold off

subplot(4,1,3,'Position',[0.15 0.3 0.8 0.2])
hold on
scatter(x2,y2,3,'filled','MarkerFaceColor',[0.4 0 0.6])
plot(x2,y2,'Linewidth',1,'Color',[0.4 0 0.6])
xlim([-13,13])
ylim([-3,3])
xticks(-12:4:12)
xticklabels([])
yticks(-2:1:2)
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
%set(gca,'XMinorTick','on','YMinorTick','on')
dummyh = line(nan, nan, 'Linestyle', 'none', 'Marker', 'none', 'Color', 'none')
legend(dummyh, '40 K')
legend boxoff
hold off

subplot(4,1,2,'Position',[0.15 0.5 0.8 0.2])
hold on
scatter(x3,y3,3,'filled','MarkerFaceColor',[0.6 0 0.4])
plot(x3,y3,'Linewidth',1,'Color',[0.6 0 0.4])
xlim([-13,13])
ylim([-3,3])
xticks(-12:4:12)
xticklabels([])
yticks(-2:1:2)
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
%set(gca,'XMinorTick','on','YMinorTick','on')
dummyh = line(nan, nan, 'Linestyle', 'none', 'Marker', 'none', 'Color', 'none')
legend(dummyh, '50 K')
legend boxoff
hold off

subplot(4,1,1,'Position',[0.15 0.7 0.8 0.2])
hold on
scatter(x4,y4,3,'filled','MarkerFaceColor',[0.8 0 0.2])
plot(x4,y4,'Linewidth',1,'Color',[0.8 0 0.2])
xlim([-13,13])
ylim([-3,3])
xticks(-12:4:12)
xticklabels([])
yticks(-2:1:2)
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
%set(gca,'XMinorTick','on','YMinorTick','on')
dummyh = line(nan, nan, 'Linestyle', 'none', 'Marker', 'none', 'Color', 'none')
legend(dummyh, '60 K')
legend boxoff
hold off

set(gcf,'Position',[0, 0, 200, 435])